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1. Introduction 

While the vacuum properties of most hadrons are known to reasonable accuracy nowadays, it is a 
heavily debated question how these properties change inside nuclear matter. In particular, for the 
light vector mesons, various theoretical predictions regarding their in-medium properties have been 
suggested. 

Among these expected in-medium effects, a so-called "collisional broadening" of the meson spec- 
tral function, due to inelastic collisions with the hadronic medium, is presumed. A second class of 
predictions claims that the vector meson masses will be shifted in the medium. These changes of 
the peak mass are connected to the partial restoration of chiral symmetry in the medium and have 
been studied via QCD sum rules. This effect has been claimed to be seen in experiments, but is 
still being discussed controversially. For a recent review on in-medium effects, see [1]. 
For studying in-medium effects, the more prominent hadronic decay modes of the vector mesons 
are unfavorable, since they are affected by strong final-state interactions with the hadronic medium 
- in contrast to the rare dilepton decay modes, which only feel the electromagnetic force. Therefore 
the latter are ideally suited to carry the in-medium information outside to the detector, undisturbed 
by the hadronic medium. 

Experimentally, dilepton spectra from elementary nuclear reactions are being studied for example 
with the CLAS detector at JLAB, where photons with energies of a few GeV interact with nuclei 
[2], or by the E325 experiment at KEK, where 12 GeV protons are used as projectile [3]. Also, 
the HADES detector at GSI has an ambitious program for measuring dilepton spectra from p+p, 
p+A and A+A reactions [4]. On the side of the hadronic decays, most notably CO — > 71 y is being 
investigated by the CB/TAPS group in photon-induced reactions at the ELS A accelerator [5]. 
In this paper, we apply the GieBen Boltzmann-Uehling-Uhlenbeck transport model (GiBUU) [6] 
to the p+p and p+Nb reactions studied by the HADES collaboration. We use GiBUU to gener- 
ate dilepton events and pass them through the HADES acceptance filter, in order to compare our 
calculations directly to the experimental data measured by HADES. 

2. The GiBUU Transport Model 

Our tool for the numerical simulation of dilepton spectra is the GiBUU hadronic transport model, 
which provides a unified framework for various types of elementary reactions on nuclei as well 
as heavy-ion collisions [6] . This model takes care of the correct transport-theoretical description 
of the hadronic degrees of freedom in nuclear reactions, including the propagation, elastic and 
inelastic collisions and decays of particles. The GiBUU model is based on the Boltzmann-Uehling- 
Uhlenbeck equation, 

(d 1 +f p H-V r -V r H-V p )f i (7,p,t)=I co n[f i ,f j ,...} , (2.1) 

which describes the space-time evolution of the one -particle phase-space density /; of a given 
particle species i under the influence of a mean-field potential. The phase-space densities /; at time 
t depend on the spatial coordinates r and the four momentum p. The right-hand side of the equation 
is given by the collision term 7 co n, which describes collisions and decays of particles. The left-hand 
side, the so-called "Vlasov" part, describes the propagation of particles in a mean field, where H is 
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the relativistic one-particle Hamiltonian. In order to solve the BUU equation numerically, we rely 
on the test-particle ansatz. Here the phase-space densities are approximated by a large number of 
test particles, each represented by a 8 -distribution in coordinate and momentum space. 
The propagation of particles with density-dependent spectral functions (usually referred to as "off- 
shell propagation") poses a particular challenge. Our approach to this problem is based on the off- 
shell equations of motion of test particles, as given in [7]. Such an off-shell treatment is necessary 
for including in-medium modifications of the spectral functions (e.g. collisional broadening of the 
vector mesons). The collisional width inside a nuclear medium of density p can be related to the 
collision cross section Onx with the nucleons in low-density approximation as 



where v re i is the relative velocity and the brackets indicate an integration over the Fermi momentum 
of the nucleons. In general this collisional width will depend on the momentum of the involved 
particle X. However, a consistent treatment of such a momentum-dependent collisional width is not 
possible with the off-shell equations of motion mentioned above, as it will lead to superluminous 
test particles. Therefore we have to neglect the momentum dependence and keep only the linear 
density dependence 



where po = 0.168fm is the normal nuclear matter density. The value of To should on average 
match the momentum-dependent width as obtained from the collision term. For the p and ft) 
mesons we typically use To = 150MeV. 

3. The Collision Term / Elementary Cross Sections 

The collision term contains all sorts of scattering and decay processes: elastic and inelastic two- 
body collisions, decays of unstable resonances and even three -body collisions (which are only 
relevant at high densities). 

The two-body part of the collision term is separated into two different regimes in terms of the 
available energy y/s : a resonance model description at low energies and the Pythia string model 
at high energies. 

For baryon-baryon collisions, the transition between the two is usually performed at y/s" = 2.6GeV. 
There is a small window around this border (±0.2 GeV), where both models are faded linearly into 
each other in order to ensure a smooth transition. For meson-baryon collisions, the transition region 
lies at V? = 2.2±0.2GeV. 

The low-energy part is given by a resonance model [8], where basically all collision cross sections 
are assumed to be dominated by the excitation of nucleon resonances. The GiBUU model currently 
contains 61 baryon species (around 30 nucleon resonances plus strange and charmed baryons). The 
resonance model gives a good description of exclusive channels like one- and two-pion production 
as well as single 77 and p production. 

At higher energies the resonance model breaks down. There we rely on the Monte Carlo event 
generator Pythia (v6.4.24) [9, 10], which is based on the Lund string model. Although Pythia 



r c0 U = P (VrelOiVx) 



(2.2) 



r co ii — To 

Po 



(2.3) 
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clearly has its strengths in the high-energy regime, we use it down to energies of a few GeV. This 
works surprisingly well, as has recently been demonstrated for example by GiBUU's successful 
description of pion data measured by the HARP collaboration [11]. In the following, we will show 
that our approach also provides a rather good description of HADES dilepton spectra. 
Since GiBUU is intended to be a multi-purpose generator, we follow the philosophy that Pythia 
is used with most of its default settings. It is obvious, however, that tuning a few selected param- 
eters can lead to a much better agreement with specific experimental data; cf. appendix A for the 
parameter set used for the HADES calculations. 

Furthermore, we had to extend PYTHIA in some more aspects. For instance PYTHIA 's treatment of 
spectral functions is rather simplified and not sufficient for a proper description of dilepton spec- 
tra, which are very sensitive to the vector meson spectral functions. Therefore we had to replace 
Pythia's internal spectral functions by proper ones, because the former lack mass-dependent de- 
cay widths and are restricted to a narrow region around the pole mass. This also gives us the 
opportunity to introduce density-dependent spectral functions. 

4. Dilepton Decays and Form Factors 

In the GiBUU model the following dilepton decay modes are taken into account: 

• direct decays, as V — > e + e~ with V = p°, CO, <f> or T] — > e + e~ , 

• Dalitz decays, as P — >■ e + e~ y with P = n Q , 7], tj' or CO — > n°e + e~ or A— >Ne + e~ . 

Most of them are treated similarly as in [12]. The leptonic decay widths of the vector mesons are 
taken under the assumption of strict vector-meson dominance (VMD), 

4 

rw e -Gu) = c V7i f (4.i) 

with the constants Cy listed in table 1 (taken from [13]). 



V 


ray(MeV) 


r ee (keV) 


Cy = T ee /my 


p 


775.49 


7.04 


9.078 • 10~ 6 


CO 


782.65 


0.60 


7.666 -10- 7 





1019.455 


1.27 


1.246 -10" 6 



Table 1: Dilepton-decay constants for V — > e + e . 

While the direct decay of the r\ meson into a i u + /i _ pair has been observed, for the corresponding 
e + e~ decay only an upper limit of BR(tj — > e + e~) < 2.7 • 10~ 5 is known [14]. However, the 
theoretical expectation from helicity suppression is still four orders of magnitude lower [15]. 
The Dalitz decays of the pseudoscalar mesons, P = 71°, 77 , T/ are treated via the parametrization 
[16], 
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with r>^ yr = 7.8 • l(T 6 MeV, T^yy = 4.6 • l(T 4 MeV and the form factors 

Fjfl(n) = 1 +b K oLL 2 , b n o = 5.5 GeV- 2 , (4.3) 

F?j(iU)= ( 1_ Af) ' A n=°- 676GeV • ( 4 - 4 ) 

The above value of has been recently determined from the HADES data at 2.2 GeV beam 
energy [17]. It should be noted that the form factors of the n° and rj Dalitz decays are sufficiently 
constrained by data, while the experimental constraints of the v\' form factor are much weaker 
[16]. A VMD form factor for the 7]' Dalitz decay can be found for example in [18]. However, 
the rj' contribution to the HADES dilepton spectra turns out to be practically insignificant. The 
parametrization of the ft) Dalitz decay 



dr. 



CO— viVe^e 



l^(M)| 2 



2oc r^^^oy 

3k jit 

At 



1 + 



-mi 



(a4 



in 



2\2 



3/2 



(4.5) 



( A 2-^)2 +A 2 r 2 



0.65 GeV and r 



is adopted from [19, 12] with r ffl ^ 3r o r = 0.703 MeV, A m 
here that the form factor of the ft) Dalitz decay is also well-constrained by data [20]. 
For the A-Dalitz decay, we use the parametrization from [21], 



(4.6) 

75 MeV. We note 



dr^Ne+e- 2a 

dpi ~3^i r ^ Nr > (4 - 7) 

r = « K+m/v) 2 [{mA+mN)2 _^ ] l/2 t(mA _ m;v) 2_ M 2 ] 3/2 |FA(M)|2) (4g) 

16 m^ntff 

where we neglect the electron mass. The electromagnetic N-A transition form factor Fa(a0 i s an 
issue of ongoing debate. Unlike the other semileptonic Dalitz decays, it is poorly constrained by 
data. At least at the real-photon point {}X = 0) it is fixed by the decay A — > Ny to |Fa(0)| = 3.029, 
and also in the space-like region this form factor is well-constrained by electron scattering data 
on the nucleon. However, this form factor is basically unknown in the time-like regime, which is 
being probed by the A Dalitz decay. 

From the theoretical side, many parametrizations are available for the space-like part, but most of 
them are not applicable in the time-like region. One of the few models which take care of the 
continuation to the time-like region is the two-component quark model given in [22] . 



5. Dilepton Spectra from Elementary p+p Collisions 

It is very important to make sure that one understands the elementary reactions before moving on 
to heavier systems, which involve effects of the nuclear medium. Fortunately, HADES has also 
measured dilepton spectra from elementary p+p reactions. These provide a base line for exploring 
the heavier nuclear systems. 
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Figure 1: (Color online) Dilepton spectra for p+p@3.5GeV. a) Mass spectrum without A form factor, b) 
mass spectrum with A form factor. 

Fig. 1 shows a comparison plot of a GiBUU simulation to HADES data [23] for a proton beam 
of 3.5 GeV kinetic energy impinging on a fixed proton target. This setup corresponds to a center- 
of-mass energy of = 3. 18 GeV. The theoretical results have been corrected for the HADES 
acceptance and reasonably reproduce the shape of the data over most of the mass spectrum. In the 
intermediate mass region around 500-600 MeV it seems like the inclusion of a proper transition 
form factor for the A Dalitz decay [22] (shaded area) is crucial for describing the data. Without 
such a form factor, the calculation strongly underestimates the experimental data in this region (by 
at least a factor of two). Another channel which could possibly contribute in the intermediate mass 
region is 77 — > e + e~. The current upper limit for the branching ratio of this decay would overshoot 
the HADES data by at least a factor of four. However, there is no r\ peak visible in the data, 
and also the theoretical expectations from helicity suppression are still orders of magnitude below 
the current experimental limit. Therefore it is rather unlikely that the direct 77 decay would give 
significant contributions to the HADES dilepton spectrum at 3.5 GeV. 

It should be noted that the default Pythia parameters already give decent 71° and X] production 
cross sections, while the vector meson production is overestimated quite a bit. This is cured by 
using the tuned parameters from appendix A. 

In order to understand the underlying processes, it is not sufficient to consider only the mass spec- 
trum. Other observables can give further insight into the reaction dynamics and can serve as a cross 
check for the validation of theoretical models. We choose to examine the transverse momentum 
and rapidity distribution in three different mass bins (see fig. 2): 

1. The low mass bin (m < 150MeV) is clearly dominated by the 7t° Dalitz channel, with some 
admixture of the r\ Dalitz. 

2. The intermediate mass range of 150 MeV < m < 550 MeV mostly contains contributions of 
the 77 and A Dalitz decays. 

3. The high mass bin (550 MeV < m) is where the vector mesons are situated. 
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Figure 2: (Color online) a) Transverse momentum and b) rapidity spectra of dilepton pairs for three mass 
bins (with A form factor). 

Distinguishing several mass bins is useful in order to separate the contributions of different chan- 
nels. In all the three mass bins, we achieve a rather good agreement with preliminary data (not 
shown) [24], in transverse momentum as well as rapidity. With Pythia's default parameters, the 
Pt spectra show a slight deviation from the data. The rapidity spectra turned out to be rather in- 
sensitive to the details of the reaction dynamics and to most of the parameters we modified. Their 
shape is mostly governed by the geometrical acceptance of the detector, which we apply to our 
simulations via the HADES acceptance filter (HAFT, version 2.0) [25, 26]. Therefore, the agree- 
ment in the rapidity spectra is mostly useful as a consistency check for the filtering process. The 
overall agreement proves to be quite good, although there are minor deviations at forward rapidity. 
The asymmetry of the deviation indicates that the deviation is indeed caused by acceptance filtering 
problems, since the unfiltered physical distributions are symmetric around mid-rapidity, and any 
asymmetries can only be introduced by the acceptance filter. 

Within this level of agreement in the elementary p+p collisions, we have a good baseline for study- 
ing in-medium effects in p+Nb, although the issue of the A Dalitz form factor is not completely 
settled. 

6. Dilepton Spectra from p+Nb Collisions 

In p+Nb reactions there will be couple of additional effects, compared to the elementary p+p re- 
actions. First of all, the primary p+N collisions will be nearly identical, apart from binding effects 
and some Fermi smearing, but besides p+p also p+n collisions will play a role. Furthermore, the 
produced particles will undergo final state interaction within the Nb nucleus, and processes like 
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Figure 3: (Color online) Dilepton mass spectra for p+Nb@3.5GeV. Left column: without A form factor. 
Right column: with A form factor. From top to bottom: vacuum spectral functions for the vector mesons, 
collisional broadening, 16% mass shift, collisional broadening plus mass shift. 



meson absorption and regeneration may become important. The secondary collisions will at aver- 
age have lower energies than the primary N+N collisions, therefore also the low-energy resonance 
part of the collision term will be involved. Finally, also the vector meson spectral functions may be 
modified in the nuclear medium. 

Fig. 3 shows simulated dilepton spectra for p+Nb collisions at 3.5 GeV in various scenarios, com- 
pared to the data from [27]. The overall agreement is not quite as good as in the p+p case. Already 
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Figure 4: (Color online) pj and rapidity spectra for three mass bins: m < 150 MeV, 150 MeV < m < 
550 MeV and 550 MeV < m. 

in the pion channel we slightly overestimate the data. This might have various reasons: too little 
absorption or too much secondary pion production in GiBUU, or even a normalization problem in 
the data. 

According to [28], the data have been normalized by comparing charged pion spectra measured 
by HADES in p+Nb to those measured by the HARP collaboration. However, the cross sections 
obtained by HARP had to be extrapolated to the slightly different beam energy and nuclear target of 
HADES. This procedure is responsible for most of the systematic error of the data (roughly 28%), 
which is shown as a gray band in the figures. 

Another striking feature of the p+Nb system is that the simulation gets close to the data in the 
intermediate mass range, even without any A form factor. Including the form factor will slightly 
overshoot the data. It seems that most of the intermediate-mass gap observed in p+p is filled up 
by low-mass p mesons, which presumably are produced in secondary collisions. Even in p+p 
collisions, one might already get a similar effect by describing p meson production via resonance 
excitation, which could give stronger contributions in the low-mass part of the p spectral function 
than Pythia's string fragmentation model. 

The mass spectrum above 500 MeV can receive further modifications from the inclusion of in- 
medium effects in the vector meson spectral functions. Fig. 3 shows a few typical in-medium 
scenarios: The first one includes a collisionally broadened in-medium width as described in chapter 
2, while the second one assumes a pole mass shift according to 

m*(p)=moU-a£\ , (6.1) 



9 



In-Medium Properties of Vector Mesons in a Transport Approach 



Janus Weil 



with a scaling parameter a = 16%. The third scenario combines both of these effects. The modifi- 
cations introduced by these scenarios are roughly on the same order of magnitude as the systematic 
errors of the data. This fact, together with the discrepancy in the n° channel and the uncertainty of 
the A form factor, presently disables us to draw any hard conclusion on vector meson properties in 
cold nuclear matter from the HADES data. 

The shape of the pj and rapidity spectra, depicted in fig. 4, show a similarly good agreement to 
preliminary data (not shown) [29] as in the p+p case. The shown pj and rapidity spectra do neither 
include a A form factor nor any in-medium effects for the vector mesons and are not significantly 
sensitive to such modifications. 

7. Conclusions 

We have shown that the HADES data from elementary p+p collisions at 3.5 GeV kinetic beam 
energy can be described very well by the Pythia event generator with a few adjusted parameters, 
as employed by the GiBUU model. It was found that the intermediate mass region can be explained 
by including a VMD-like transition form factor for the A Dalitz decay. Also the p+Nb reaction at 
the same beam energy is reasonably well described by the GiBUU transport model. However, a 
few discrepancies are left, which prevent us from drawing final conclusions on in-medium effects 
at this point. We infer that the HADES data at 3.5 GeV can not be understood without fixing the 
open issue of the A Dalitz form factor. 

A. Pythia Parameters 

GiBUU relies on the Pythia event generator for elementary collisions at cm. energies above 2.2 
(2.6) GeV for meson-baryon (baryon-baryon) systems. Pythia's default parameters are tuned to 
rather high energies, just like most of the available parameter sets (called 'tunes' in PYTHIA). To 
apply Pythia successfully to collisions at HADES energies, we adjusted a couple of parameters. 
The following parameter set represents an attempt to create a 'HADES tune' for Pythia. A 
thorough description of all the parameters can be found in the Pythia manual [10]. 



Parameter 


Default 


Professor 


Our 




value 


tune [30] 


value 


PARJ(ll) 


0.5 


0.31 


0.15 


PARJ(12) 


0.6 


0.4 


0.2 


PARJ(21) 


0.36 


0.313 


0.25 


PARJ(25) 


1.0 


0.63 


0.63 


PARJ(26) 


0.4 


0.12 


0.12 


PARP(91) 


2.0 




0.25 



Table 2: Pythia parameters tuned to HADES energies. 

The most prominent modification is the suppression of the vector meson contributions by tuning 
down PARJ(1 1) and PARJ(12) by a significant amount. PARJ(25) and PARJ(26) are the parameters 



10 



In-Medium Properties of Vector Mesons in a Transport Approach 



Janus Weil 



for suppressing the 77 and v\' contribution, which we adopt from the 'Professor' tune [30]. More- 
over, the amount of transverse momentum has to be adjusted via PARP(91) and PARJ(21), which 
control the primordial kj and the pj from fragmentation, respectively. 
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